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Abstract 

^\j \ We study the binding energies of singlet trions, i.e. charged excitons, in carbon nanotubes. 

The problem is modeled, through the effective-mass model, as a three-particle complex on the 
surface of a cylinder, which we investigate using both one- and two-dimensional expansions of the 
wave function. The effects of dimensionality and correlation are studied in detail. We find that 



qh' the Hartree-Fock approximation significantly underestimates the trion binding energy. Combined 

with band structures calculated using a non-orthogonal nearest neighbour tight binding model, the 
^' , results from the cylinder model are used to compute physical binding energies for a wide selection 
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of carbon nanotubes. In addition, the dependence on dielectric screening is examined. Our findings 
indicate that trions are detectable at room temperature in carbon nanotubes with radius below 



I. INTRODUCTION 

Many applications of semiconducting carbon nanotubes (CNTs) rely on the formation 
of junctions using either dopedi or asymmetrically gated CNTs^. Such junctions are the 
basis for e.g. CNT based light emitting diodes^''^'^ and other optoelectronic devices. In the 
description of the optical properties of CNTs, excitons are known to significantly influence 
emission^'^ and absorption^. The importance of excitons in CNTs is a consequence of the 
extremely large binding energy of several hundred meV-"-. Thus, in an optically excited 
junction with injected carriers, an exciton in combination with either a hole or an electron 
may form a three-particle complex known as a trion. It has been shown experimentally 
that trions form in biased GaAs quantum wells^ under optical excitation via the interaction 
between excitons and injected electrons and holes. Consequently, similar processes are ex- 
pected to occur in CNTs. Hence, if the binding energy is sufficiently large, these complexes 
are expected to be of importance for optoelectronic devices and could even be detectable at 
room temperature. 

In traditional semiconductors, trion binding energies and optical spectra have been stud- 
ied in one^.'iS, two^i"— and three dimensionsi^. In contrast, trions in CNTs have only been 
studied in an effectively one- dimensional modeUS, in which carriers were assumed to be com- 
pletely delocalised around the circumference of the tube. Within this model, it was found 
that singlet and triplet trions in CNTs are stable against dissociation into excitons and free 
carriers. Furthermore, the difference in energy between the positive 5^ and negative S~ 
singlet trion states was rather small for mass fractions a = rrie/mh G [0.8; 1], where me is 
the effective electron mass and rrih is the effective hole mass. Also, numerical trion binding 
energies were estimated for CNTs with radius from 2.5A to lOA. It is clear, however, that 
a one-dimensional model may have limited applicability for large diameter CNTs, in which 
motion around the circumference is an important degree of freedom. 



The one-dimensional model used in Ref. [ig has many potential applications as it greatly 



simplifies the complexity of the calculations for trions, biexcitons, charged biexcitons, etc. 
However, it still remains to assess the quality of this approximation and, therefore, it is 
important to compare one- and two-dimensional solutions for trions. It is clear that the 



effective-mass model used in Ref. Il6| is an approximation to more accurate ab initio ap- 
proaches. However, comparison with ab initio results for the exciton energy and its scaling 
with tube radius^ as well as the exciton wave function^^ demonstrates excellent agreement 
between these different approaches^^ii^. In this paper, we investigate two simplifications 
of the full effective-mass trion problem: First, the one-dimensional model is examined by 
comparison to two-dimensional solutions of the full problem. Second, the accuracy of the 
Hartree-Fock approximation for trion binding energies is determined in both one- and two- 
dimensional models. As our starting point for both models, we expand the wave function 
in a convenient basis. The same basis expansion is applied in the Hartree-Fock calculation, 
which is briefly outlined. Next, we investigate the binding energies of trions, in both the 
one- and two-dimensional models, as functions of cylinder radius r and mass fraction a. We 
compare the Hartree-Fock solutions with the full solutions (at a = 0.0), and we analyze the 
distribution of electrons along the circumference in a negative trion state S~ for the two 
methods. Finally, the CNT band structure is used to convert model results into physical 
binding energies for a wide range of CNTs, using both the one- and two-dimensional models. 

II. THEORY 

Adopting the fundamental trion equation from Ref. il5i and introducing relative cylindrical 
coordinates along with removal of center-of-mass motion, we get the Hamiltonian for the 
negative trion 

where r is the nanotube radius and Xi and rOi denote axial and circumference coordinates, 
respectively, as shown in Fig. [1] The exciton operators, hi and /i2, are given by 

and the Coulomb potential in the cylinder geometry is given as^ 

v{x,e) = ^= ^ (3) 
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Here, the coordinate pairs (xi,6'i) and (x2,^2) describe the motion of the first and second 
electron relative to the hole, respectively. The positive trion Hamiltonian Hj^ is found by 




FIG. 1. Geometry of a negative trion on the surface of a cylinder of radius r, described in relative 
coordinates. The filled circles illustrate electrons and the open circle illustrates a hole. 

the replacement a — t- a^^. Note that H^ and i^+ have incorrectly been interchanged in Ref. 
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The Hamiltonian is expressed in natural exciton units, i.e. effective Bohr radii a*^ and 



effective Rydbergs Ry* for distances and energies, respectively. 



If the radius r is much smaller than the effective Bohr radius a^, the electrons are highly 
delocalised around the circumference and the angular dependence of the wave function for 
the lowest states will be nearly constant. In contrast, when r is comparable to a^, the 
angular dependence is expected to have a cusp-like behaviour similar to the one found for 
excitons^. A function that suitably describes this behaviour is the absolute value of sine 
I sin||. Thus, it would be reasonable to model the angular part of the wave function as a 
linear combination of the two limits. Moreover, to model the electron-electron repulsion for 
H^ and hole-hole repulsion for if+, the function | sin( ^^^^^ )| constitutes a sensible choice as 
it attains its minimum along the diagonal 9\ = 92- Having this in mind, we expand the wave 



function according to 

Ni,N2,N3 4 

ij{Xi,9i,X2,92) = ^ ^Cijkl(l)i{Xi)(f)j{x2)(f)k{Xl~ X2)(pi{0i,92), (4) 

i,j,k I 

with Gaussians 4>i{x) = e~"^^ along the translational direction and an angular basis {(j)i}f^i 
along the circumference, where 
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Here, i, j and k run from 1 to A^^i, A^2 and N^, respectively, and N = Ni ■ N2 ■ N3 is the total 
number of Gaussians used in the expansion. The expansion coefficients Cijki are obtained 
from the eigenvector c of a matrix equation of the form {K + U)-c = E^ S ■ c. We define 
the trion binding energy Eb = Ex — Ex as the difference between the exciton binding energy 
Ex and the trion energy E^, where a positive Eb indicates a stable trion. The kinetic matrix 

'^ijkl,i'j'k'l' 



K and overlap matrix S elements turn out as Sijki,i'j'k'i' = SiJ^,,,yS\i, and 
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where 5ij is the Kronecker delta and with the angular overlap integrals S^, given in Tab. 
m In the derivation, care should be taken as the second derivative of the absolute value 
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TABLE I. Table of angular overlap integrals. 

of a sine gives a delta function, which takes many of the kinetic energy elements to zero. 
The potential energy matrix elements turn out to be slightly more complicated and can be 
expressed in terms of 4 different Meijer G functions. 

In order to formulate the problem in terms of the Hartree-Fock approximation we now assume 
that the eigenfunction of the operator Eq. ([1]) is a Slater determinant of single electron states 
x{x, 0) with anti-symmetry in the spin part. From standard quantum mechanics^^ we then 
obtain the Fock operator for the corresponding system 
^ a^ 1 92 2 



Qrj,2 ^2 QQ2 



+ Vh{x 
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where the Hartree potential Vh is given by 
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From the eigenvalue equation Fx{x, 9) = ex{x, 0) the Hartree-Fock eigenvalues are found 
and the lowest one Sq is related to the trion ground state energy in the usual way 



E\ 



2£n- 



dx / de\x{x,e)YVH{x,e). 



(9) 



Expanding x(x, 9) in a basis of Gaussians along the translational direction and absolute value 
of sines along the circumference, as done for the full problem, we are capable of solving the 
Hartree-Fock eigenvalue equation self-consistently. 

III. RESULTS AND DISCUSSION 



Optimisation of the Gaussian coefficients was carried out using a steepest descent 



method. The following coefficients, in units of a*^ , were found for the different opti- 
misations: ai e {0.143,1.16,4.98,29.0,250} for both the ID and 2D exciton and ai,aj,ak G 
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FIG. 2. Comparison between the one- and two-dimensional models. For each of the two models 
three cases are shown, namely, S^g, S-^^ q and S^g. The positive and negative trion binding energies 
were calculated for a = 0.80. The S^ q binding energies were calculated for a = 1.0. 

{0.0651,0.145,1.68,9.65,48.7} for the ID trion. For the 2D trion, the Gaussian coefficients 
were found to be ai,aj G {0.165,1.68,9.65,48.7} for the functions e"""*^ and e~"^^ , and 



ak e {0.0000171,1.68,9.98,48.7} for the functions e 



-Okix-yf 



Similarly, the coefficients for the 



Hartree-Fock calculations were found to be a^, aj E {0.0648, 0.195, 1.04, 5.28, 27.5, 99.3, 250} 
for both the ID and 2D cases. All optimisations were carried out at r^ = O.la^ and all 
coefficients were subsequently scaled with r^/r"^ for calculations involving other radii. 



The binding energies Eb of trions were found as functions of radius for both the one- 
and two-dimensional models. In Fig. [2], it is seen that the one- and two-dimensional results 
deviate less as r approaches zero, as expected. It is important to realize that the effective 
Bohr radius a*^ is dependent on nanotube species via the effective masses and dielectric 
constant e. Due to the scaling of masses with nanotube radius, it turns out that r is, in fact, 
roughly proportional to a^. For CNTs embedded in a medium with e = 3.5 it is found that^ 
r ~ O.la^, and for this value the correction is roughly 13% of the two-dimensional result. 
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FIG. 3. Trion binding energy as a function of the mass fraction a = me/rrih for the one- and 
two-dimensional models. The calculations were done for r = O.la^. 



The one- dimensional model performs significantly worse when the radius is increased. For 
r ^ 0.3a^ the one-dimensional result is 42% below the two-dimensional one and, as would 
be expected, the one-dimensional model can only be considered a good approximation for r 
close to zero. Reahstically, though, r > 0.2a^ will only occur for CNTs embedded in media 
with small screening e < 1.75. 



From Eq. ([T]) it follows that the mass fraction a affects the binding energy via a mixed 
kinetic energy term coupling the two relative coordinates (xi,^i) and (a;2;^2)- The mass 
fraction a of CNTs varies between 0.86 and 1.0^^ and to gain an idea of the upper and 
lower bounds of the trion binding energies we have plotted the positive and negative trions 
at cr = 0.8 in Fig. [2l Actual binding energies should lie between the Sq^ and 5*^0 energy 
curves and the S^^^ and Sq^^ curves for the negative and positive trion, respectively. It is 
seen that the difference is fairly small for both the one- and two-dimensional models. In Fig. 
131 we illustrate the binding energies as functions of the mass fraction a for both models for 
r = O.la^. In both cases, the binding energy for the negative trion S~ was found to have 
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FIG. 4. Solutions for tiie one- and two-dimensional models compared with the corresponding 
Hartree-Fock solutions. The calculations were carried out at a = 0.0. 

little variation for a G [0; 1]. As a consequence, the results calculated for S^q are, in fact, 
excellent approximations for both S'^ and S~ with mass fractions in the range a G [0.8; 1]. 
It follows that neglecting the mixed kinetic energy term is an acceptable approximation, 
which in the worst case would give an error of 3.2% at r = O.la^. 



Next we calculated the binding energies using the Hartree-Fock approximation. As 
demonstrated above, the mixed kinetic energy term is only a minor correction for actual 
nanotubes. Thus, in the following we ignore this term and limit the discussion to Sqq states. 
In Fig. m the Hartree-Fock energies are plotted along with the Sqq trion energies for both 
one- and two-dimensional solutions. The Hartree-Fock solutions significantly underestimate 
trion binding energies is both cases. Differences up to a factor of 2, in the two-dimensional 
case, were found. The best Hartree-Fock result obtained in the interval r G [0; 0.3] was at 
the end point r = 0.3a^ where the Hartree-Fock energy is 60% of the correct result. The 
trion problem of two electrons interacting with a positive hole is mathematically similar to 
the He atom. Hence, the dramatic failure of the Hartree-Fock approximation is surprising 



at first sight. However, the failure is readily ascribed to several characteristic factors in the 
trion case. Primarily, we focus here on the binding energy Eb = Ex — Et given as the 
difference between exciton and trion energies. Hence, relatively small errors in the Hartree- 
Fock trion energies Ej^ lead to large errors in the binding energies provided Ej, — E^ is 
comparable to Ex — Et- Secondly, Fig. H] shows that the error increases as we decrease the 
radius, i.e. when the effective dimension of the system decreases from two towards unity. 
Actual CNTs can be regarded as effectively 1.7 dimensional systems^^. The approximate 
treatment of electron-electron repulsion in the Hartree-Fock calculation is more serious in 
a low-dimensional geometry due to the increased overlap between electrons enforced by the 
geometry. Hence, the Hartree-Fock error for the (three dimensional) He atom is small even 
if the energy is measured relative to the hydrogen atom. Finally, in the present problem, 
the strengths (charges) of the electron-hole and electron-electron are identical. In He, the 
larger nuclear charge Z = 2 enhances the electron- nucleus interaction. As a result, the 
contribution to the energy from the electron-electron repulsion is relatively smaller and, 
thus, the Hartree-Fock approximation is less inaccurate. 

To investigate the behaviour of the wave function around the circumference the proba- 
bility distribution, averaged over the translational coordinates xi and X2, as a function of 
angles, 

/•OO /"OO 

Pt(^1,^2)= / / dXidX2\Mxi,ei,X2,92)? 
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AAA 



= 22z2z2 ^ijk,i'j'k'(^ijklCi'j'k>l'Vl{di,62)Vv{0i,92), 
1,1' i,j,k i',j',k' 

has been shown in the left column of Fig. [5^-c for r = 0.05a^, r = O.lOa^ and r = 0.25a^. 
The probability distribution gives the probability PT{di,62)d6id62 of finding the electrons 
with angles between 9i and 9i + d9i, and, 02 and 02 + d92 from the hole, respectively. It is 
normalised such that integration over both angles from — vr to vr yields 1. It is seen that the 
probability distribution tends to be more delocalised as r tends to zero. This is also what we 
would expect since a purely one-dimensional system has no angular dependence. Further it 
is noticed that the distribution has a butterfly shape as is seen in other casesF^'^. In Fig. [Sjl 
we have plotted the exciton probability distribution as function of the angle for r between 
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0.05a^ and r = 0.25a^. The exciton probability distribution is given by 

/oo -^ 5 j 

dx\^Pxix,e)\' = $^$^Q,c,,, /^^^f )(^)^W(e), 
-oo • ., • •, \ Oil -\- Oil' 

where ip\ (6) = l/-\/2vrr and y^g (9) = \ sin ||/v^27rr. This plot confirms the contracting 
behaviour of the trion probability distribution seen in a-c. In Fig. |5fe, the Hartree-Fock prob- 
ability distribution has been plotted for r = O.lOa^. Clearly, the effects of approximating 
the electron-electron repulsion are largest along the diagonal 9i = 02- This is emphasised in 
the difference plot Fig. [5]F, where the difference between the Hartree-Fock distribution and 
the full distribution has been plotted in percent. The errors in the Hartree-Fock distribution 
vary from -2% up to 4%. The delocalization of the trion wave function is enhanced by the 
nearly equal electron and hole masses. Hence, the centre of mass will not coincide with the 
position of any of the constituents but, rather, lie somewhere in between. This is another 
notable difference between trions and He atoms, for which the centre of mass will practically 
fall on top of the nucleus. As a consequence, both electrons and holes are almost completely 
delocalized around the circumference in the realistic trion case a ~ 1.0 and r ^ O.la^, as 
illustrated in Fig. [5]d. 

Whether trions in CNTs will be detectable at room temperature is determined by the 
magnitude of the trion binding energy Eb relative to the thermal energy kBT. The trion 
binding energy depends on the static dielectric constant e as well as the reduced electron- 
hole pair mass fi, through the effective Rydberg Ry* = 13.6eV ■ /i/e^ and the effective 
Bohr radius a*^ = 0.529A ■ e/^^^'^^. The reduced mass /i is obtained from the CNT band 
structure, and in order to determine the hole and electron masses in CNTs, for a particular 
chiral index {n,m), we used a non-orthogonal nearest neighbour tight binding model. The 
transfer integral was chosen as^^ t = — 2.89eV and the overlap as s = 0.1. The soundness 
of these parameters follows from the fact that the predicted Fermi velocity of graphene 
vp = 9.6 ■ lO^m/s is within 5.7% of the average experimental value^^. Since the CNT band 
structures are derived from the graphene band structure this indicates that effective hole 
and electron masses will be estimated with a similar error. Using this tight binding model 
the reduced mass was found as /i = l/{m^^ + m~^). The dielectric constant e on the other 
hand, is determined partly by the CNTs and partly by their surroundings. It has been shown 
in several articles^^""— that the contribution from the surroundings plays a significant role in 
the determination of optical transition energy in CNTs. Further, the electro-static potential 
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FIG. 5. Probability distributions as a functions of angles for (a) r = 0.05a^, (b) r = O.lOa^ and (c) 
r = 0.25a^ for the singlet trion Si q calculated with a basis expansion. The exciton probability 
distributions have been shown in (d) for the radius varying from r = 0.05a^ to r = 0.25a^. In 
(e), the Hartree-Fock probability distribution is shown at r = O.la^ and (f) shows percent wise 
difference between the Hartree-Fock solution and the solution obtained from the full problem at 
r = 0.1a|j 
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0, in an quasi one-dimensional nano structure confined along y and z, will be governed by 
the perpendicular component e±, since^^ oc {ez^yX^ + exSzy'^ + ^x£yz'^)~^^'^ ~ ^I^l^l""*^ if the 
tube radius is small. It has been shown that e^, is small for a wide range of nanotubes in 
Ref. I29I, and therefore the surrounding material is expected to dominate the static dielectric 
constant e. Moreover, the results expressed in effective units are universal results and do 
not depend on e. The dependence arises only through the conversion into physical units (eV 
and A) via the effective Rydberg Ry* and Bohr radius a^, similarly to the excitoi>i^ and 
biexciton^. Since the dielectric constant depends on the material in which the CNTs are 
embedded, a*^ and Ry* should be determined for the individual experiment. Therefore e can 
be regarded as an experimental parameter varying from sample to sample. Well-established 
methods for synthesising samples with high concentrations of (6,5) CNTs exist^>2^ and we 
will therefore use (6,5) as an example in the following. CNTs having chiral index (6,5) 
have in some cases been suspended in sodium cholate^ and, in other cases, embedded in a 
polymeric matrix^^. The dielectric constant e should be chosen according to either of these 
materials. However, no exact values of e were found for the above suspension materials. It 
has been shown that the photonic transition energy as a function of environmental dielectric 
constant saturates near— Senv ~ 5.0. Expecting that the surrounding material will have 
a dielectric constant somewhat above the dielectric constant of air, e = 3.5 is considered 
as a reasonable average. The radius for (6,5) CNTs is r = 3.73A and the effective masses 
of electron and hole turn out to be rrie = 0.0803 and rrih = 0.0866, respectively, which 
leads to a reduced mass /x = 0.0417. The effective Rydberg and Bohr radius are found to 
be Ry* = 0.0462eV and a*^ = 44. 5A, respectively. With these values the effective radius 
is r = 0.084a^, and using this, the negative trion binding energy can be calculated to 
Eb = 1.28Ry* ~ 59meV. Hence, for this species Eb is larger than fc^T = 26meV and the 
trion state for (6,5) CNTs is expected to be detectable at room temperature. 



In the hope that these results will stimulate measurements of the trion binding, we have 
calculated the energies at e G [2.0; 5.0], which corresponds to finding the binding energy for 
(6,5) CNTs suspended in a large variety of solutions. We found that the binding energy, in 
this range for e, will be at least Eb = 36meV and at most Eb = 132meV. The result has 
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FIG. 6. Binding energy of negative trions in (6,5) CNTs as a function of dielectric constant e. 

been shown in Fig. [6] in a double logarithmic plot and we also include a curve fit given by 

EB{e) ^ (0.372£"^-5*^ + 0.00608)meV. (10) 

The result is not surprising as it simply emphasises the importance of choice of suspen- 
sion material with respect to measuring the binding energy of trions. As a consequence, 
experiments should be conducted with suspension materials with low dielectric constants. 
Also, the power dependence is what would be expected as a similar result was found for 



the exciton^, where it was shown that Ex^ 



r oc r 



-0.6 



Ry*, with r in units of a*^. Using 



the effective Rydberg and Bohr radius it is easily shown that the exciton energy for CNTs 
follows a similar power law Ex{s) oc 6~'^'^. 

Keeping e = 3.5 constant, we found the positive and negative trion binding energies for 
all semiconducting CNTs with SA < r < ISA. The results are seen in Fig. [71 The binding 
energies were found using both the one- and two-dimensional models. As seen in the inset, 
the two-dimensional model has improved the energy, compared with the one-dimensional 
one, with up to 15%, and on average 11%. It is noticed that in order to observe trions in 
CNTs at room temperature, the radius of the CNTs must be below approximately SA. Since 
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FIG. 7. Binding energy of trions in CNTs in physical units using effective electron and hole masses 
derived from a non-orthogonal tight binding scheme. The shaded area illustrates the instability 
region for which Eb < ksT. Inset: percent-wise difference between the one- and two-dimensional 
models. 

the binding energy increases with decreasing radius, CNTs with low radius constitute better 
candidates for observing trions. Hence, CNTs with chiral index (6,5) are very promising 
candidates for measuring trion binding energies. 



IV. CONCLUSION 

In this work, trions in CNTs have been modeled as three-particle complexes bound to 
the surface of a cylinder. We have shown that the angular behaviour plays a significant role 
with a contribution of 13% to the binding energy for cylinders with radius r = O.lOa^. It is 
concluded that the energy for the Sqq trion constitutes a fairly good approximation to both 
positive and negative singlet trions with mass fraction above 0.80. We have demonstrated 
that the Hartree-Fock method applied to the Sqq trion equation yields results that at best, 
are 60% of the correct result. We conclude that the binding energies of trions in CNTs are 
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lowered with 11% on average by including the angular part of the wave function for trions. 
Finally, trions are expected to be detectable in doped CNTs with r < 8A and we consider 
CNTs with chiral index (6, 5) as a very good candidate for measuring trion binding energies. 
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